TWO  SPACE-CONSTRAINED  MULTI-ITEM 
DETERMINISTIC  INVENTORY  MODELS. 


Klaus  Lieding 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  tCMOOL 

rONTEREY.CAUFORWAWMO 


QTP 


it  Ml 


P" 


Monterey,  California 


TH 


TWO  SPACE-CONSTRAINED  MULTI-ITEM 
DETERMINISTIC  INVENTORY  MODELS 

by 

Klaus  Lieding 

September  1975 


Thesis  Advisor: 


A.W.    McMasters 


Approved  for  public  release;  distribution  unlimited. 


TH9785 


UNCLASSIFIED 


SECURITY   CLASSIFICATION  OF    THIS  PAGE  (When  Date  1 

Entered) 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.     REPORT   NUMBER 

2.  GOVT   ACCESSION  NO. 

3.     RECIPIENT'S  CAT  ALOG  NUMBER 

4.     TITLE  (and  Subtlllm) 

Two   Space-Constrained   Multi-Item 
Deterministic    Inventory   Models 

5.     TYPE  OF   REPORT  ft   PERIOD  COVERED 

Master's    Thesis; 
September    1975 

«.     PERFORMING  ORG.    REPORT   NUMBER 

7.    authorc*; 

Klaus    Lieding 

8.     CONTRACT  OR  GRANT  NUMBER^ 

9.     PERFORMING  ORGANIZATION   NAME   ANO  ADDRESS 

Naval   Postgraduate   School 
Monterey,    California      93940 

10.     PROGRAM   ELEMENT.  PROJECT,   TASK 
AREA  ft   WORK  UNIT  NUMBERS 

11.     CONTROLLING  OFFICE  NAME   AND   ADDRESS 

Naval   Postgraduate   School 
Monterey,    California      93940 

12.     REPORT   DATE 

September    1975 

13.     NUMBER  OF   PAGES 

60 

14.     MONITORING  AGENCY   NAME  ft    ADDRESSf//  dlllerent  from  Controlling  Office) 

IS.     SECURITY  CLASS,  (ol  thle  report) 

Unclassified 

15».     DECLASSIFI  CATION/ DOWN  GRADING 
SCHEDULE 

16.     DISTRIBUTION   STATEMENT  (ol  thlt  Report) 

Approved   for   public   release;    distribution   unlimited. 

17.     DISTRIBUTION  STATEMENT  (ot  the  abatrmct  entered  In  Block  20,  It  dlllerent  from  Report) 

16.     SUPPLEMENTARY  NOTES 

19.     KEY  WORDS  (Continue  on  reverae  elde  II  neceeeary  and  Identity  by  block  number) 

20.     ABSTRACT  (Continue  on  reverae  elde  II  neceeeary  and  Identity  by  block  number) 

This   paper   presents    two   types   of   multi-item  deterministic 
constrained   inventory  models   which   relate   to   the    constraints 
of    the   equivalent   of   a   direct    support   unit   in   the   Federal 
German   Army;    namely,    that   there    is    only    a    limited   amount   of 
space   available    for    storage   of    items.      One    is   a   continuous 
review  model  with   no   stockouts    allowed   and   the   other   is    a 
periodic    review  model  where   backorders    are   allowed.      In    contrast 

FORM         ii?1 

1  JAN  71    1473 


DD 

(Page    1) 


EDITION  OF    1  NOV  88  IS  OBSOLETE 
S/N    0102-014- 6601   | 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGe  (When  Data  Entered) 


UNCLASSIFIED 


JtCUt^TY    CL  ASSIFtC  »  TlON    OF    Thi',  P»r;E'H>#n    n,im   F.nl»r*3 


(20.   ABSTRACT   Continued) 

to  previous  models  analyzed  by  Naddor  and  others,  time 
phasing  of  the  arrivals  of  orders  of  different  items  is  allowed 
The  benefit  derived  from  this  phasing  is  the  reduction  in  the 
average  amounts  of  on-hand  inventory  and  backorders .   Solutions 
are  obtained  for  the  two-item  case  and  the  n-item  case  is 
discussed.   Examples  illustrate  the  results. 


DD      Form       1473 

,  1  Jan  73 
S/N    0102-014-0601 


UNCLASSIFIED 


SECuniTv  classification  of  this  PAGErr***"  £>•'•  emtrmd) 


Two  Space-Constrained  Multi-Item 
Deterministic  Inventory  Models 

by 

Klaus  Lieding 
Captain,  Federal  German  Army 
Betriebswirt  (grad) ,  FHS  I  of  the  Army,  1970 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September  1975 


ABSTRACT 

This  paper  presents  two  types  of  multi-item  deterministic 
constrained  inventory  models  which  relate  to  the  constraints 
of  the  equivalent  of  a  direct  support  unit  in  the  Federal 
German  Army;  namely,  that  there  is  only  a  limited  amount  of 
space  available  for  storage  of  items.   One  is  a  continuous 
review  model  with  no  stockouts  allowed  and  the  other  is  a 
periodic  review  model  where  backorders  are  allowed.   In 
contrast  to  previous  models  analyzed  by  Naddor  and  others, 
time  phasing  of  the  arrivals  of  orders  of  different  items  is 
allowed.   The  benefit  derived  from  this  phasing  is  the 
reduction  in  the  average  amounts  of  on-hand  inventory  and 
backorders.   Solutions  are  obtained  for  the  two-item  case 
and  the  n-item  case  is  discussed.   Examples  illustrate  the 
results . 
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I.   INTRODUCTION 

Military  forces  have  to  be  supplied  with  all  kinds  of 
items  in  order  to  complete  their  missions.   Because  many  of 
the  needed  items  cannot  be  purchased  from  outside  suppliers 
at  the  time  the  forces  need  them,  the  items  have  to  be 
purchased  and  stored  earlier  than  actually  needed. 

From  the  point  of  view  of  the  ultimate  user  it  would  be 
best  to  store  all  such  needed  goods  as  near  as  possible  to 
his  own  location,  so  that  whenever  an  item  is  needed  it  can 
be  delivered  immediately.   However,  there  are  obviously 
constraints  which  make  this  course  of  action  impossible. 

First,  there  is  a  budget  constraint  which  allows  only  a 
certain  amount  of  money  to  be  spent  on  purchases  of  equipment, 
spare  parts  and  other  needed  goods  and,  second,  there  is  a 
space  constraint  because  the  supply  support  units  are  required 
to  be  mobile  so  that  they  can  stay  in  close  proximity  to  the 
combat  units  they  support.   In  the  case  of  direct  support 
units  for  land  forces  the  budget  constraint  is  either  non- 
existent or  overwhelmed  by  the  space  constraint. 

In  the  Federal  German  Army  these  units  are  called  BVP's 
and  each  consits  of  a  fleet  of  trucks.   They  are  equivalent 
to  the  Direct  Support  Units  (DSU's)  of  the  U.  S.  Army.   The 
BVP's  constitute  the  last  echelon  of  the  German  Army's  supply 
system  and  are  not  involved  in  buying  goods  from  outside 
suppliers.   Therefore  the  budget  constraint  on  purchases  is 


nonexistent.   However,  the  number  of  available  trucks  for 
a  BVP  is  limited  by  its  operational  requirements  and,  to  a 
certain  extent,  by  the  available  budget  for  the  armed 
services.   Therefore  any  relevant  inventory  model  has  to  take 
a  space  constraint  into  consideration  when  calculating  optimal 
operating  policies. 

Two  models  will  be  analyzed  in  this  thesis.   The  first 
will  be  a  multi-item  deterministic  lot  size  model  subject  to 
a  space  restriction  which  requires  that  the  amount  in 
inventory  at  any  time  cannot  exceed  the  available  storage 
space  and  a  restriction  that  the  time  between  orders  for  all 
items  must  be  the  same.   The  variables  which  are  being 
optimized  in  this  system  are  the  time  between  orders  for  each 
item  and  the  times  between  orders  for  different  items.   The 
notion  of  allowing  orders  for  different  items  to  arrive  at 
different  times  will  be  referred  to  as  time  phasing  of  the 
orders  for  different  items.   The  solution  for  this  model  is 
obtained  for  the  two-item  case  and  the  procedure  for  solving 
the  n-item  case  is  outlined.   The  second  model  is  a  multi-item 
deterministic  order  level  system.   It  is  subject  to  the  same 
space  and  cycle  time  constraints  as  the  lot  size  model. 
However,  stockouts  are  allowed.   The  decision  variables  are 
the  order  levels  of  each  of  the  items  and  the  time  phasing 
of  the  ordering  of  each  item.   This  model  is  solved  for  the 
two-item  case  and  discussed  for  the  n-item  case. 


II.   PREVIOUS  WORK  IN  THE  AREA 

A.   CONSTRAINED  LOT  SIZE  MODEL 

A  multi-item  lot  size  model  with  a  single  space  constraint 
has  been  examined  by  Churchman  [2]  and  Holt  [6] .   The  model 
was  defined  in  the  usual  form  of  deterministic  inventory 
problems;  that  is,  to  minimize  the  total  variable  costs  — 
reorder  plus  carrying  costs  —  of  the  system.   The  space 
constraint  required  that  the  sum  of  the  average  space 
requirements  for  all  items  not  exceed  the  total  space  available 
for  these  items.   The  assumptions  which  were  made  in  the 
formulation  of  the  problem  were: 

1.  a  known  constant  demand  rate; 

2.  infinite  replenishing  rate; 

3.  constant  lead  time; 

4.  shortages  are  not  allowed. 

The  cycle  lengths  of  all  items  were  not  required  to  be 
identical.   In  addition,  the  time  phasing  of  each  item  was 
not  addressed. 

A  disadvantage  of  this  model  is  that  only  the  sum  of  the 
average  inventories  must  be  no  greater  than  the  space  available. 
This  requirement  can  be  satisfied  mathematically  and  yet  not 
be  met  in  reality  because  the  possible  different  cycle  lengths 
may  cause  the  order  for  more  than  one  item  to  arrive  at  the 
same  time. 


B.   A  CONSTRAINED  ORDER  LEVEL  SYSTEM 

A  multi-item  deterministic  order  level  model  has  been 
examined  by  Naddor  [7] .   This  model  used  the  first  three 
assumptions  of  the  above  model  by  Holt  and,  in  addition, 
required  a  fixed  review  interval  for  all  items.   Backorders 
were  allowed  however.   The  objective  function  for  minimization 
was  the  total  variable  cost  —  carrying  plus  backorder  costs  — 
and  the  constraint  required  that  the  summation  of  the  space 
requirements  of  the  optimum  order  levels  not  exceed  the  total 
space  available. 

The  major  difference  between  Naddor 's  model  and  the  one 
to  be  presented  in  this  paper  is  that  in  his  model  all  orders 
arrive  at  the  same  time,  while  the  model  being  developed  here 
permits  time  phasing  of  arrivals  of  orders  of  different  items. 
His  solution  for  the  optimal  order  level  for  the  n-item  case 
is: 


c 


S,  =  q 


n  .  -  g  b. 
2i    ^o  1 


1    ^i  c,  .  +  c~ . 
li     2i 


where 


S.  =  order  level  of  item  i  when  lead  time  is  zero? 

q.  =  replenishment  size  of  item  i; 
c, .  =  the  carrying  cost  of  item  i  per  unit  of  time; 
Cj-    ~   the  shortage  cost  of  item  i  per  unit  of  time; 

b.  =  unit  space  required  for  item  i; 
B  =  total  warehouse  space  available; 
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and 


n 
=      (    I 


c~ .b.q. 
2i    1^1 


1=1      li  2i 


-    B) 


n 
Z 


b.2q. 

1       ^1 


-    c,  .    +   c„ . 
1=1      li  2i 
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III.   A  CONSTRAINED  MULTI-ITEM  LOT  SIZE  SYSTEM 

The  model  which  is  developed  in  this  chapter  is  an  n  item 
lot  size  model  constrained  by  the  amount  of  space  which  is 
available  for  storage  of  inventory.   The  following  assumptions 
are  made : 

1.  Demand  for  item  i  is  deterministic  at  a  constant 
positive  rate  of  r.  units  per  unit  time. 

2.  No  shortages  are  allowed  and  the  stock  of  an  item  will 
be  immediately  replenished  whenever  its  inventory  level 
reaches  zero. 

3.  The  lead  time  is  constant. 

4.  The  replenishment  rate  is  infinite. 

5.  One  type  of  storage  is  usable  by  all  items. 

6.  The  inventory  carrying  cost  is  proportional  to  the 
average  amount  in  inventory  at  any  time;  its  dimension 
is  $  per  unit-time. 

7.  The  order  cost  is  a  constant  which  is  incurred  each 
time  an  order  is  placed;  its  dimension  is  $. 

The  variables  which  will  be  optimized  to  produce  minimum 

cost  are  the  length  T  of  the  reorder  period  for  every  item 

and  the  time  x.  at  which  an  order  for  the  ith  item  arrives 

after  arrival  of  an  order  of  item  number  1  (x-,  =  0)  .   The  cost 

equation,  as  a  function  of  the  scheduling  period  for  any 

single  item  under  these  assumptions,  is  given  by 


Ci(T)  =  cliriT/2  +  C3i/T  ;  (1) 


where     c.  .  =  the  inventory  carrying  cost,  and 

c->.  =  the  unit  reorder  cost. 
3i 
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This  notation  is  that  of  Ref.  7,  page  58.   For  n  items,  the 
total  cost  can  be  expressed  as 


n 

TC   =    I   (c  .r.T/2  +  c,./T)   , 
i=l    X1  x        J1 


given  that  all  items  are  using  the  same  reorder  period  length 
The  space  constraint  can  be  formulated  as 


n 

Z   b.I.(t)  <  B      for  all  t  >  0  ,  (2) 

i=l   x  1 


where 

B     =  total  space  available; 

b.    =  the  space  occupied  by  a  unit  of  the  ith  item; 

1^ (t)=  the  inventory  level  of  the  ith  item  at  time  t. 

Since  maximum  inventory  levels  occur  at  the  times  at  which 
orders  are  received,  these  times  correspond  to  critical  points 
relative  to  the  constraint.   Therefore  the  constraint  needs 
only  to  be  examined  at  these  arrival  times.   This  will  result 
in  n  constraints  for  an  n-item  case  because  of  the  possibility 
that  each  item  may  arrive  at  a  different  time. 

To  solve  the  problem  let  A.  represent  the  ith  Lagrange 
multiplier.   Then  the  Lagrangian  becomes 


n  n 

L   =    I      C.  (T)   +    £   X.F.  (3) 

1=1  1=1 
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where  C. (T)  is  given  by  equation  (1)  and  F.  represents  the 
ith  constraint  function  (describing  (2)  at  each  x. , 
i=l,2, . . . ,n) .   The  solution  to  equation  (3)  will  be  presented 
for  cases  which  have  two  and  n  items.   The  two-item  case  is 
presented  to  illustrate  the  general  solution  procedure. 

A.   SOLUTION  OF  THE  TWO-ITEM  CASE 

Equation  (3)  for  the  two-item  case  reduces  to: 


L   =   c11r1T/2  +  c31/T  +c12r2T/2  +  c32/T  +  ^-^l   +  A2F2  * 


Typical  fluctuations  of  the  inventory  levels  of  the  items  can 
be  represented  as  shown  in  Figure  1. 

From  Figure  2  it  can  be  seen  that  the  constraint  function 
at  time  zero  when  item  number  1  arrives  will  be 


b,r,T  +  b2r2x2  <_  B   , 


and  that  the  constraint  function  at  time  x2  when  item  number 
2  arrives  will  be 


b1r,(T  -  x2)  +  b2r2T  <_  B 


Conversion  of  inequalities  into  equalities  for  use  in  the 
Lagrangian  yields: 

Fl  "  blrlT  +  b2r2X2  +  Vl2  "  B   ;  (4) 

F2  =  t^r^T  -  x2)  +  b2r2T  +  V^   -   B   ;  (5) 

where  V,  and  V2  are  slack  variables. 


r](T-x2) 


*  t 


FIGURE  1:     Fluctuation  of  inventory  levels  for  two  items  for 
the  deterministic  lot  size  model. 
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I(t)'k 


b-jr-jT  =  b2r2T 


b^TpXp  =  b,r,  (T-Xp) 


FIGURE  2:     Space  occupied  by  two  items  when     b^r-j   =  b2r\. 


16 


Combining  (3),  (4)  and  (5)  gives 


L   =   clir]T/2  +  c31/T  +  c12r2T/2  +  c^/H 


+  A^b-j^T  +  b2r2x2  +  V^   -  B) 

+  X2(b1r1(T-x2)  +  b2r2T  +  V22  -  B)  .   (6) 


Taking  partial  derivatives  of  the  Lagrangian  with  respect  to 
V1  ,  V"2  /  x2  /  T,  X,  and  X2,  setting  them  equal  to  zero,  and 
simplifying,  yields: 


X2y1  =   0       ;  (7) 


X2V2  =   0  (8) 


Alb2r2  "  A2blrl      =   °       '  (9) 


Cllrl/2  "  C31/T2  +  C12r2/2  "  C32/t2  +  Xlblrl 


+  A2blrl  +  ^2b2r2   =   0   ;    (10) 


blrlT  +  b2r2x2  +  Vl      "  B   =   °  '*  (13L) 


b1r1(T-x2)  +  b2r2T  +  V ^   -   B  =   0      .  (12) 


Equations  (7),  (8)  and  (9)  indicate  that  there  are  only 
two  possible  cases: 
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1.  A,  jt   0   and   A2  ^  0  which  implies  that  V,  =  V2  =  0; 

2 .  A^  =  A2  =  0   ,   and  the  V^  values  can  be  arbitrary. 

When  A,  ^   0  and  A2  ^  0 ,  both  constraints  are  binding.   Solving 
equations  (11)  and  (12)  with  V,  =  V"2  =  0  for  T  and  x~ 
respectively,  yields: 


T*    =  B(b1r1 +b2r2)/(b12r12  +b1r1b2r2   +   b22r22)     .  (13) 


*  2      2  2      2 

x2    =  Bb2r2/(b1   r,       +  t,,r,b2r2    +  b2    r2    )     .  (14) 


Substituting   T   and   x2    into   the    total    cost   equation  yields 
the  minimum  cost  C     where 


Co    =  B(cllrl  +  C12r2)  (blrl  +b2r2)/2D  +    D(C31  +  C32)/A    ' 


and 


D  =  b12r12  +  birib2r2  +  b^r/   ; 


A  =  B(b1r1  +  b2r2) 


The  ratio 


x*         b2r2 


blrl  +  b2r2 
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provides  a  convenient  means  of  comparing  x~  to  T  for  this 
completely  constrained  case.   If  b,r,  =  b2r.?  then 


*  * 

4  =  \    .   If   blrx  <  b2r2  then   4  >  \   .   If  b^  >  b2r2 

* 

X2    1 
then  — 5f-  <  y  . 

T    * 


For  the  unconstrained  case,  equation  (9)  requires  that 
if  one  of  the  X   values  equals  zero,  the  other  has  to  be  zero 
too.   The  optimal  unconstrained  value  of  T,  denoted  as  T  , 
can  be  determined  from  equation  (10) .   It  is 


2(c   +  c  2) 
T   =   I  ££__ £i .  (15) 

1  Cllrl  +  C12r2 


Equation  (15)  is  the  same  expression  as  derived  by  Naddor  [7]. 
The  value  for  x~  is  dependent,  however,  on  the  value  of  B. 
The  value  of  B  which  corresponds  to  the  dividing  point  between 
cases  1  and  2  is  obtained  by  equating  equations  (13)  and  (15) 
and  solving  for  B.   This  value  is  denoted  as  B,,  and  is  given 
by  equation  (16) 


-1  =  V^J  Tu      •  <"» 


When  B  _<  B,  ,  equations  (13)  and  (14)  give   the  values  of  T' 

and  x^  .   When  B  >  B,  ,  T   =  T   and  the  value  of  x0  is  no 
2  1  '        u  2 

longer  restricted  to  that  given  by  (14)  . 
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If  B  >  B2  where  B2  =  (b1r1+b2r2)  Tu  '   then  x2  can  be 
any  value  in  the  interval  0  <_  x2  <  T  .   This  corresponds  to 
B  being  so  large  that  even  if  orders  for  both  items  arrive 
simultaneously  they  will  not  be  faced  with  a  storage 
restriction.   When  B,  <  B  <  B2  ,  the  storage  constraint  will 
not  be  restrictive  provided  that  x2  is  so  selected  that  the 
constraint  is  not  exceeded  at  the  time  either  item  arrives. 

When  the  first  item  arrives  (at  t  =  0)  the  total  combined 

on-hand  inventory  consumes   b,r,T  +  b~r~x~   amount  of 

J  1  1  u    2  2  2 

storage  space.   When  the  second  item  arrives  (at  t  =  x_)  the 

total  combined  on-hand  inventory  consumes  space  in  the  amount 

of  b,r,  (T  -  x_)  +  b~r0T   .   The  value  of  x~  must  be  selected 
1  1   u    2      22u  2 

so  that  the  total  combined  on-hand  inventory  at  either  time 
does  not  exceed  B.   Thus  if  x2  lies  in  the  interval 


B~  -  B  B  -B 

-T7E-±  x2  1  Tu  "  TF    '  (17) 

rlDl     ^     U     r2  2 


then  both  conditions  will  be  satisfied  when  B,  <  B  <  B2 . 
When  B  =  B„  the  interval  given  by  (17)  simplifies  to 
0  —   x2  —  Tu  *   As  B  decreases  below  B2 ,    the  upper  bound 
decreases  and  the  lower  bound  increases  at  the  same  rate. 
The  interval  reduces  to  a  single  value  for  x2  at  B  =  B,; 
the  corresponding  x?  value  is  given  by  equation  (14).   As  a 
consequence  of  this  behavior  of  the  interval  for  B  between 
B,  and  B2 ,  it  follows  that  equation  (14)  gives  an  optimal 
value  of  x2  for  all  B  >  B,  ,   and  the  only  optimal  value  for 
B  <_  B1. 
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B.   SOLUTION  OF  THE  n-ITEM  CASE 

The  Lagrangian  for  n  items  is: 


n  n 

L  =   E   (c  .r.T/2  +  c  ./T)  +   £   A  F.  ,  (18) 

i=l    ±1  1        J1      i=l   x  1 


with  the  constraint  functions  given  below.   These  functions 

assume   0  <  x~  <  x->  <  •  •  •  <  x 
—  2  —  3  —     —  n 


blrlT  +  b2r2*2  +  b3r3x3  +         •-         +bnrn^  +  Vl2  =  B 


b^fl^)  +  b2r2T  +  b3r3(x3-x2)  +   ...   +  Vn(V*2J  +  V2*  =  B 


Vl^3  +  t>2r2(T-  (x3-X2))  +  b3r3T  +  ...  +  b^x^)  +  V*     =  B 


bLrl(%-l,+V2(T-(Vfx2,)+     •••+bnrn(VVl)+Vn-2l     =     B     . 


blrl^V  +  b2r2 (T  -  ^"x2)  +       '  *  *        +  bnrnT  +  Vn*  "  B 


After  forming  the  Lagrangian,  taking  partial  derivatives  with 
respect  to  x2  ,  x..,  ...,  x   and  V,,  V"2  ,  ...,  V  ,  and  setting 
them  equal  to  zero,  it  can  be  shown  that  either  all  of  the 
Lagrange  multipliers  are  zero  or  all  slack  variables  are  zero, 
The  partial  derivatives  with  respect  to  A,,  A2 ,  ...,  A  , 
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when  the  V. 's  are  all  zero  give  back  equality  constraints  at 
t  =  0,x9,...,x   .   Solving  these  equations  for  T,  x0,    .  ..,  x  , 
yields  : 


n 
B  l      b.r. 
i=l   x  x 


(19) 


B 


l 

Z 

k=2 


bkrk 


X      1-  f  <^  f   *  •   •   f  ll      f 


(20) 


where 


n    9   ?    n-1         n 
D  =   z   b.  r.   +   Z   b.r.    z   b.r.  . 
i=l   1   x    i=l   x  L  j=i+l   D  J 


(21) 


The  resulting  value  of  C  ,  the  minimum  total  variable  cost, 
is : 


n 


n 


C0  =  B  (.Z1  Cliri/2  (blrl+'"Vn)/D  +  D  <*  c31^  (blrl+"-4bnrn). 
i=l  1=1 


In  the  unconstrained  case,  optimal  T  is  obtained  from 


u 


.1 

1 

2      Z      c 
i=l      J1 

! 
i 

| 

n 

.E,  curi 

1=1 

(22) 


The  n-item  equivalents  of  B,  and  B2  are: 


D  T 


B. 


u 


n 

Z   b.r. 
i=l   x  x 


(23) 
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B2    "    .".    biriTu  •  <24> 

1  =  1 


And,  if  B  _>  B„ ,  then  T  =  T   and  each  x.  can  take  on  any 

value  in  the  interval  0  <_   x.  <_  T   .   When  B,  <  B  <  B2  then 

* 

the  x. 's  must  be  chosen  in  a  way  that  the  constraints  are 

not  violated.   As  B  decreases  below  B~  ,  the  upper  bound 
decreases  and  the  lower  bound  increases.   However,  equation 
(20)  has  been  obtained  for  the  case  where  the  constraints 
have  to  be  met,  i.e.  where  we  had  equalities  in  all  of  the 
constraint  functions.   Changing  the  equalities  to  inequalities 
will  result  in  an  infinite  set  of  optimal  feasible  solutions. 
But  x.  as  given  by  equation  (20)  will  always  be  in  this  set 
of  optimal  feasible  solutions  because  the  solution  to  the 
more  restricted  problem  will  always  be  an  answer  to  the  less 
restricted  problem. 

The  set  of  optimal  solutions  creates  the  ranges  for  the 
x.  values.  They  could  be  obtained  using  linear  programming 
in  the  following  way: 

1.  Possible  objective  functions  are 

a)  min  x. ,  i  =  2,3,...,n,  to  get  the  smallest  lower 
bound  for  x. ; 

b)  max  x.,  i  =  2,3,...,n,  to  get  the  greatest  upper 
bound  for  x. . 

2.  The  constraints  are  all  of  those  shown  below 
equation  (18) . 
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For  solving  problems  which  have  n  >  2  items  the  flow  chart 
of  Figure  3  is  helpful  and  summarizes  the  preceding  analysis. 

C.   EXAMPLES 

Two  examples  are  presented  below  to  illustrate  the  model. 

Example  1 

Consider  a  warehouse  which  stores  two  different  items 
whose  parameters  are: 

Item  1:  Item  2 : 

r,   =  200  units  per  year  r2  =  250  units  per  year 

b1      =    5  ft3/unit  b   =    3  ft3/unit 

c-ji  =  $  2  per  unit-year  ci2=  ^  ^  Per  uni^Y631" 

c31  =  $  30  c32=  $  25 

If  total  available  warehouse  space  is  400  cubic  feet,  what 

* 

are  the  optimal  order  period  T   for  each  item  and  the  optimal 

arrival  time  x2  for  item  number  2  after  an  order  of  item 
number  1  arrives? 

Calculating  T   and  Bn  yields: 

3   u      1  J 

Tu  =  0.41  years   ;     B1  =  543  ft3   . 

Because  B  <  B,  ,  both  storage  constraints  are  binding. 

*  * 

Therefore  we  get  T  =0.30  years  and  x2  =  0.13  years  by 

using  equations  (13)  and  (14). 
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calculate: 

T  =  equation  (22) 

B-,  =  equation  (23) 


yes 


I  no 

calculate: 
B?  =  equation  (24) 

T*  =  equation  (19) 

* 
x.  =  equation  (20) 


1 

no 
f 

* 

T     =  T 
u 

* 
x,   =  equation  (20) 

FIGURE  3:  Flowchart  for  solving  lot  size  model  problems 
having  n  >  2  items. 
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Example  2 

Suppose  the  data  of  example  number  1  remains  unchanged 
except  for  B.   Assume  B  is  now  600  ft  . 

The  values  for  Tu  and  B-j^  are  unchanged.   But  since  B  is  now 
greater  than  B,,  evaluation  of  B2  is  required.   It  is 

calculated  to  be  719  ft3. 

* 

Because  B1    <  B  <  B2  ,   T   =  T   =0.41  years  and  x2  lies 

in  the  interval  given  by  (17);  that  is, 
0.12  years  <_   x2  <_     0.25  years  . 

Checking  the  feasibility  at  time  t  =  0  shows  that,  using  the 

* 
lower  bound  as  the  x2  value,  the  occupied  warehouse  space 

3  * 

would  be  500  ft  ;  using  the  upper  bound  as  the  x2  value, 

3 

600  ft  would  be  needed.   The  respective  space  needs  for 

*  3  3 

time  t  =  x?  would  be  600  ft   and  467  ft  . 

D .   COMMENTS 

The  constrained  multi-item  lot  size  system  developed  in 
this  chapter  differs  from  that  of  Holt  [6]  in  the  following 
ways  : 

1.  Holt's  model  allows  a  different  cycle  length  for  each 
item,  while  a  common  cycle  length  for  all  products  is 
used  in  the  above  model. 

2.  Holt's  constraint  required  that  the  sum  of  average 
inventories  be  no  greater  than  the  space  available,  while 
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the  constraint  function  used  in  the  above  model  requires 

that  the  space  occupied  by  the  inventory  at  any  time  be 

no  greater  than  the  total  space  available. 
3.   The  time  phasing  between  orders  for  different  items  is 

calculated  in  the  above  model ,    while  the  time  at  which 

an  order  arrives  is  not  regulated  by  Holt. 
The  second  aspect  is  perhaps  the  most  important  difference 
because  the  actual  space  restriction  in  Holt's  model  is  a 
nebulous  notion.   Realistically,  the  space  constraint  restricts 
the  maximum  on-hand  inventory  rather  than  the  average  level. 

It  is  interesting  to  note  that  the  solution  for  the 
scheduling  period  can  also  be  used  to  calculate  a  minimum 
amount  of  storage  space  needed  given  a  specific  scheduling 
period.   The  equation  for  the  warehouse  space  in  terms  of  the 
scheduling  period  is 

n 
B   =   TD/  I      b.r. 
i=l   !  1 

where  D  is  given  by  equation  (21) . 
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IV.   A  CONSTRAINED  MULTI-ITEM  ORDER  LEVEL  SYSTEM 

A.   THE  MODEL 

The  system  which  is  analyzed  in  this  chapter  is  an  n-item 
order  level  or  periodic  review  system  which  is  constrained  by 
the  total  amount  of  space  which  is  available  for  storage  of 
the  items.   The  analysis  of  this  system  will  be  conducted 
under  the  following  assumptions: 

1.  Demand  for  item  i  is  deterministic  at  a  constant  rate 
of  r.  units  per  time. 

2.  The  period  between  replenishments  or  scheduling  period 
is  a  fixed  constant  T. 

3.  The  replenishment  quantity  of  the  ith  item  raises  the 
inventory  at  the  time  the  quantity  is  received  to  the 
order  level  S-  after  backorders  are  made  up. 

4.  The  replenishment  rate  is  infinite. 

5.  The  lead  time  is  constant. 

6.  The  same  type  of  storage  is  used  by  all  items. 

7.  The  inventory  carrying  cost  is  proportional  to  the 
average  amount  in  inventory.   Its  dimension  is  $  per 
unit  time. 

8.  The  backorder  cost  is  proportional  to  the  average  shortage 
level.   Its  dimension  is  $  per  unit  time. 

9.  The  order  costs  are  ignored  because  they  are  a  fixed  cost 
if  the  scheduling  period  is  a  prescribed  constant. 

From  these  assumptions  it  can  be  seen  that  the  lot  sizes  which 

are  being  considered  are  constant  since 


*i  =  ri  T  ' 
where  q.  represents  the  ith  item  lot  size 
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The  variables  which  are  to  be  optimized  in  this  system 
are  the  order  levels  (S-  for  the  ith  item)  and  the  points  in 
time  at  which  each  item  arrives  (x.  for  the  ith  item,  where 
xx  =  0)  . 

The  cost  for  one  scheduling  period  for  the  ith  item  as  a 
function  of  the  order  level  can  be  shown  to  be  [7] : 


c(si) 


CliSi2  ,  C2i<ri-Si)2 

2r.T        2r.T       ' 


where   c, .  =  inventory  carrying  cost, 

and     c~  .  =  inventory  shortage  cost  (using  Naddor1 s 
notation) . 


Then  the  total  cost  equation  for  an  n-item  model  can  be 
expressed  as: 

2  2 

n    c.  . S  .     c0  .  (r. -S  •  ) 

CMS   S       S  )  =   E   (-Ai_i-  +   2l   x  X      ) 

1  =  1        1  1 

The  storage  constraint  can  be  formulated  as: 


n 

Z   b.I. (t)  <  B 
i=l   x  1 


for  all  t  =  0  , 


where 


B  =  total  space  available, 

).  =  amount  of  space  which  is  occupied  by  one  unit 
1   of  item  i, 


I±(t) 


inventory  level  of  the  ith  item  at  any  point 
in  time  t. 
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As  in  the  preceding  model,  the  maximum  amount  in  inventory 
of  any  item  occurs  immediately  after  the  ordered  item  arrives. 
Therefore  the  arrival  time  is  a  critical  time  with  regard  to 
the  space  constraint.   Because  a  constraint  will  be  required 
for  the  time  of  arrival  of  each  item,  there  will  be  n 
constraints  for  n  items. 

Figure  4  illustrates  the  two  possibilities  in  fluctuations 
of  inventory  for  two  items.   Two  cases  are  possible  because, 
at  the  end  of  each  scheduling  period,  an  order  for  one  item 
arrives  and  the  other  item  is  either  on  backorder  or  in  stock. 
Care  must  be  taken  in  formulating  the  constraints  because  when 
an  item  is  on  backorder  its  inventory  level  is  negative  but 
the  amount  of  its  occupied  space  is  zero. 

B.   SOLUTIONS  FOR  THE  TWO- ITEM  CASE 

The  problem  for  the  two-item  case  can  be  stated  as 
follows : 


Minimize 

2  .2,2  2 

__       C11S1         c21(rlT~Sl)  C12S2         C22(r2T"S2)  ,?c-, 

TC  "     2r][T     +  2TJT  2r2T  2r^  '"  U^ 


subject   to   the   constraints 


b^   +   max[0,    b2  (S2  -  r£  (T-x2)  )  ]    <_  B         ;  (26) 


b2S2    +   max[0,    b1(S1-r1x2)]  <   B         .  (27) 
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Case  1 


Case  2 


I(t)' 

i 

sl 

&2. 

s. 

"X 

^ 

i 
i 

V^ 

^V. 

^g 

T 

FIGURE  4:     Fluctuation  of  inventory  levels   for  two  items   for 
the  deterministic  order  level  model. 
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Constraint  (26)  reflects  the  combination  of  the  constraint 
cases  at  t  =  0.   Whenever  S2  -r2(T-x2)  <  0  ,   which  means  that 
item  number  2  is  on  backorder,  this  constraint  reduces  to  the 
restriction  that  b,S,  be  less  than  or  equal  to  B.   Similarly, 
constraint  (27)  reflects  the  combination  of  cases  at  t  =  x«. 

1.   Totally  Unconstrained  Case 

If  the  problem  is  unconstrained,  then  the  optimal  value 

for  S-,  and  S2  can  be  obtained  by  simply  taking  the  partial 

derivatives  of  equation  (25)  with  respect  to  S,  and  S2,  setting 

them  equal  to  zero  and  solving  for  S,  and  S2.   The  solutions, 

denoted  as  S  ,  and  S  ~,  are: 
ul      u2 

C21rlT 

S  ,  =  -=~f= ;  (28) 

ul   c11+c21 


C22r2 

S    =  _££_f ;  (29) 

u2    c12  +  c22 


which  are  the  same  expressions  as  in  Naddor's  unconstrained  * 
case  [7].   Because  B  _>  b,S  -i  +k>2s  ?  E  B3  '   even  simultaneous 
arrivals  are  not  space   constrained  and  x2  can  assume  any 
value  in  the  interval   0  <_  x2  £  T  . 
2 .   Partially  Unconstrained  Case 

When  B  <  B-.  ,   then  simultaneous  arrivals  would 
violate  the  space  constraint.   If  B  is  not  too  restrictive, 
the  problem  could  still  be  unrestricted  relative  to  S  -  and/or 
S  2  if  careful  time  phasing  of  the  orders  is  done. 
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The  individual  constraints  were: 


b;,S1  +  max[0,  b2  (S2  -  r2  (T-x2)  )  ]  <_ 


B 


fc>2S2  +  max[0,  b1(S1~r1x2)]      <_  B 

Notice  that  S  -,  and  S  2  will  not  necessarily  violate  the 
space  constraint  if  B  >  B2  and  x2  is  chosen  carefully,  where 

B2  =  maxtb^,  b2Su2]   .  (30) 

They  will  therefore  remain  optimal. 

If  we  assume  that  both   b,S  n  <  B   and   b0S  0  <  B 

1  ul  2  u2 

then  the  functions   max[0,  b2(S  2-r2(T-x2))]   and 
max[0,  b,  (S  ■,  -  r,x2)  ]   at  t  =  0  and  t  =  x2 ,  respectively, 
could  be  positive.   At  t  =  0  the  value  of  the  max  function 
is  not  to  exceed  the  total  space  available  minus  the  space 
occupied  by  item  number  1.   A  feasible  positive  max  function 
at  t  =  0  satisfies: 

0  <   b2(Su2-r2(T-x2)J  1  B  "  blSul   ' 


or,  in  terms  of  x? , 


Sn  S0    B  -  b,S  , 

T  -  -S3.  <  x9  <  T  -  -J£  +      1    ul  (31) 

r2     2  -      r2      b2r2 
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s 

If  x0  <  T H£  then  b0(S  0-r(T-x0))  <  0  .   If  x.  is 

z  —      r~  z      uz  z        —  z 

greater  than  the  right-hand  bound  then  the  total  inventory 
at  t  =  0  would  exceed  B,  the  storage  restriction  value. 

When  t  =  x2  the  value  of  the  positive  max  function 
is  not  to  exceed  the  total  space  available  minus  the  space 
occupied  by  the  optimum  order  level  of  item  number  2.   This 
means  that 


0  <  b1(Sul-r1x2)  <  B-b2Su2  , 


or,  in  terms  of  x2 , 


S  ,    B  -  b0S  0         S  t 

-Hi £-*£  <  x0  <  -Hi  .  (32) 

rl     rlbl   "  2    rl 

s 

In  this  case  if   xn  >  ,   then   bn  (S  ,  -r,x0)  <  0 

2  —  r,  1   ul    12   — 

and  if   x2  <  left-hand  bound,  the  total  on-hand  inventory 
at  t  =  x2  would  exceed  the  storage  restriction. 

In  comparing  inequalities  (31)  and  (32)  it  is  clear 
that  the  space  constraint  will  be  only  met  if: 


S  ,    B  -  b0S  0  S  _    B  -b,S  , 

-Hi ^  <  x0  <  T  -  -Hi  +      1  ul  (33) 

rl       rlbl   -   2  -      r2      b2r2 


However,  feasibility  requires  that  the  left-hand  side  of  (33) 
be  less  than  or  equal  to  the  right-hand  side. 
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and 


The  two  bounds  of  (33)  are  equal  when  B  =  B  where 


Bx  =  B3  -  BQ   ,  (34) 


B3  =  blSul  +  b2Su2  '*  <35> 


blrlb2r2 

B0  =  [  b  r  -f-b  r   ]  T   '  (36) 

u      D±Ti   tD2r2 


Therefore,  if   B2  <_  B   <_  B  <  B3   then  the  problem  is 

* 

unconstrained  with  respect  to  S,  and  S2  provided  x-  is  chosen 

within  the  boundaries  of  inequality  (33) . 
Simplifying  (33)  results  in: 


B^  -  B  B-.  -  B 

TTlx2lT-TF       *  (37) 

Dlrl  D2r2 


It  follows  that  when  B  =  B_,  the  range  given  by  (37)  is 
0  <_   x2  <_   T  ;  when  B  <  B~  the  range  becomes  narrower  with  the 
upper  bound  decreasing  at  a  rate   l/b?r2  and  the  lower  bound 
increasing  at  a  rate   -l/b,r,   per  one  unit  decrease  in  the 
value  of  B.   We  note  in  passing  that  if   B  =  b,S  ,  =  b^S  2  , 
the  range  (33)  reduces  to 


S  ,  S  0 

jai  <  •    <  T  _  jS. 

r,  —  2  —      r~ 
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The  condition  resulting  in  B   >  max[b,S  , ,b~S  J  is  a  case 

x        1  ui   2  u2 

where  both  constraints  are  binding  for  any  B  <  B  while  the 

max  functions  are  both  still  positive  for  B  near  B  . 

x 

3 .   Constrained  Cases 

a.   B   >  B~   and   B„  <  B  <  B 
x    2         0        x 

If  Bx  >  max[b,S  1#b2S  2]   and  B  <  B   the 
constraints  when  the  max  functions  are  both  positive  are: 


b1S1  +  b2S2  -  b2r2T  +  b2r2x2  <  B   ; 


b2S2  +  b1S]_  -  bir]Lx2  <  B   . 


Subtracting  the  second  inequality  from  the  first  gives: 


*  b2r2  B0 

X2   =   [  b  r  +  b  r   ]  T   *  FT"   *  (38) 

z        Yl   2  2  11 


Thus,  x~  is  independent  of  B  provided  that  both  max  functions 

* 

remain  positive.   Substitution  of  x~  into  either  constraint 

gives : 


* 

b,S,  +  b2S2  _<  B  +  k>,r,x2    , 


_<  B  +  BQ 


Therefore  we  have  the  very  simple  relationship  that,  whatever 
values  of  S,  or  S2  are  chosen,  we  must  have 

blSl  +  b2S2  -  B  +  B0   '  (39) 
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Now  inequality  (39)  does  not  reflect  the  fact  that  the  space 

constraints  at  t  =  0  and  t  =  x2  would  be  violated  if  either 

*  * 

b,S,  >  B   or  k2S2  >  B  .   Therefore,  the  following  conditions 

must  also  be  imposed: 


blSl  -  B   ;  (40) 


b2S2  <  B   .  (41) 


The  form  of  the  resulting  constraint  region  is  illustrated 
in  Figure  5(a).   The  unconstrained  optimal  solution  (S  ,  ,S  2) 
lies  outside  this  region.   Because  the  objective  function  is 
convex  in  S,  and  S2  and  the  constraint  region  is  convex,  it 
follows  that  the  constrained  optimum  lies  on  the  boundary  of 
the  feasible  region. 

If  B  is  reduced  to  B~  we  get  the  condition  illustrated 
in  Figure  5(b)  and  when  B  <  B~  we  get  the  result  of  Figure 
5(c).   These  cases  will  be  analyzed  in  the  next  section. 

To  search  for  the  optimal  solution  when   BQ  <  B  <  B 
we  evaluate  S-,  and  S2  using  only  one  segment  of  the  boundary 
at  a  time.   We  begin  by  assuming  that  the  constrained  S,  and 
S2  lie  on  the  line  given  by  (39)  with  equality.   We  can  then 
use  (39)  to  solve  for  S2  as  a  function  of  S-.,  substitute  that 
function  into  the  objective  function,  differentiate  with 
respect  to  S,  ,  set  the  derivative  to  zero  and  solve  for  S-^. 
We  get 
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B  >B, 


b]S1   =  B 


I 

I  Feasible 
t  Region 

l 


W/ff/ll/WH/WMl 


4 


unconstrained  solution 


b,S1  +  b2S2  =  B  +  BQ 


B  =  B 


0 


mm@w&&mk 


c.     B  <  B 


0 


FIGURE  5:      Feasible  regions  for  the  two-item  deterministic 
order  level   model  where     b-,S,  +  bpS^  <  B  +  B~   . 
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*       .         c21  -  c22    [b7]    +    [-ET-]    [b7]    [s-T] 

c      _    c      -   2 2 2  u2 

S1    =   S-j^    =   ;  (42) 

C21  1     2       c22 

[^i]     +     [r-i]         [-££] 

bul     D2     bu2 


*  B+B       b    A 

S2  =  S2  =  1-57^1  "  tbj]  Sl     '  <43> 


If  either  S,  or  S2  violates  the   b-.S-.  <    B   or   k2S2  <  B 
constraints  then  the  optimal  constrained  values  lie  on  either 
of  these  latter  constraint  boundaries. 

If  S-,  is  larger  than  B/b,  then  the  largest  feasible 
S-,  value  is  B/b-,  .   This  then  will  be  the  optimal  value  of  S,  . 
Substitution  of   S,  =  B/b-,   into  ^':.  objective  function  and 
differentiating  with  respect  to  S2  will  give  S  2  when  the 
derivative  is  set  to  zero.   Therefore,  optimal  constrained 

S2  is  the  largest  value  the  constraints  will  allow  without 

* 

reducing  S,  .   This  corresponds  to  the  S2  value  at  the  inter- 

* 

section  of  (39)  and  (40);  that  is  S2  =  B0/b2 . 

* 
Similar  arguments  lead  to   S,  =  B^/b,   and 

*  * 

S2  =  B/b2   when  S„  is  greater  than  B/b2  . 

* 

In  summary,  if   BQ  <  B  <  B    then  x2  is  given  by 

*       * 
equation  (38).   To  determine  S,  and  S2 : 


1.   Calculate  S, ,  S2  . 


* 


2.   If   S,  <_   B/b,  ,   S2  <_  B/b2   then   S1  =  S^    ,      S2  =  S2  . 
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3.  If  S1    >  B/bj^  ,   then   S1  =  B/b   ,   S2  =  BQ/b1    . 

*  * 

4.  If   S2  >  B/b2  ,   then   S1  =  BQ/b1  ,   S2  =  B/b2  . 

b.   B   >  B2   and   B  _<  BQ  . 

Figure  5(b)  and  5(c)  emphasize  that  when  B  <  B0 
the  constraint  given  by  (35)  lies  outside  the  boundaries  of 
the  feasible  region  and  is  trivially  satisfied  by  any  solution 

in  the  feasible  region.   It  follows  from  the  arguments  in  the 

*  * 

preceding  section  (a)  that   S,  =  B/b,   and   S~  =  B/b2   and 

we  realize  immediately  that  both  max  functions  must  be  zero. 

The  analysis  of  the  preceding  section  showed  that 
if 


S  S 

—  <  x0  <  T -   , 

rl  "   2  -      r2 


then  both  max  functions  will  remain  zero.   Substitution  of 
*      * 


S,  and  S2  into  this  range  gives 


<  x0  <  T  -  r-^—    .  (44) 


blrl  "  2  -     b2r2 


This  range  reduces  to   x?  =  Bft/b,r,  (equation  (38))  when 

B  =  BQ ;  the  range  increases  to   0  <  x2  <  T   as  B  is  reduced 

to  zero. 

Figure  6  summarizes  the  ranges  of  optimal  x2  for 
the  case  where   B   >  B2  .   From  this  figure  it  becomes  obvious 
that   x?  =  B0/b,r,   is  optimal  regardless  of  the  value  of  B. 
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FIGURE  6:     The  range  of  optimal   x2  as  a  function  of  B 
for  the  case  where     Bx  >  B2    . 
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c.   B2  >  B    and   B  <  B2  . 


In  the  constrained  case  where   B~  >  B   ,   further 

reduction  into  cases  where   bnS  ,  >  b~S  ~  .   b~S  ~  >  b,S  , 

1  ul     2  u2      2  u2     1  ul 

and  k-jS  ,  =  b2S  2   are  necessary.   We  will  also  make  use  of 


B1    =  min[b1Sulfb2Su2]   .  (45) 

(1)   B2  =  b!sui/  B!  =  b2Su2'  and  Bl  <  B  <  B2  * 

When   B,  <  B  <  B2   then  the  first  constraint 

to  become  binding  is  that  at  t  =  0 .   Its  max  function  will  be 

zero.   An  appreciation  of  what  is  taking  place  can  be  obtained 

by  looking  at  the  x2  range  given  by  (33) .   When   B2  >  B    then 

the  upper  and  lower  bounds  do  not  reduce  to  the  same  value 

(at  B  =  B  )  before  the   B  -b,S  ,   term  in  the  upper  bound 
x  1  ul  ^^ 

goes  to  zero.   The  resulting  range  for  x2  from  (33)  is  now: 


b^S  «  S  0 

2  u2  u2  . . , . 

K  -  '  x9  <  T —  (46)  • 

blrl  ~   2  ~      r2 

* 

when   S.  =  B/b,   is  substituted  for  S  , .   The  lower  bound  is 
11  ul 

only  valid  however  as  long  as   B  >_  B,  =  b2S  2  .   The  optimal 

value  of  S.  is  S  v 
2      u2 

Inspection  of  (46)  shows  that  it  is 
independent  of  B  provided  that  B,  <  B  <  B2  .   When  the  value 
of  B  is  reduced  to  B, ,  (46)  reduces  to  (44);  that  is 


B  B 

<  x,  <  T  - 


blrl  -  2  "     b2r2   ' 
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and,  since  S2  must  be   B/fc>2   for  all  B  <  B,  ,  this  range 
must  be  valid  for  all  such  B. 

The  analysis  in  section  b  above  showed  that 
(44)  gave  feasible  x2  values  as  long  as   B  _<  Bfi  .   The  analysis 
in  this  section  then  suggests  that  B,  must  be  less  than  or 
equal  to  BQ  if  (46)  and  (44)  are  to  be  feasible  and  identical 
at  B  =  B,.   Now  the  range  of  (46)  is  feasible  if 


b~S  0        S  ^ 

<  T 


blrl  -      r2 


Collecting   like   terms   results    in 


b   r 

b0S    0    <    [    r- ±r4 ]    T 

2    u2   _        *>!*!+ *>2r2 


which  reduces  to  B,  <  B~  . 

In  summary,  if   B2  e  b,S  ,  >  B    and 

* 

B,  e  b2S  2  <  B  <  B2   then   x2   is  any  value  in  the  range  given 

by  (46)  and   S,  =  B/b,  ,   S2  =  S  2  .   When   B  <_   B,   we  again 
have   x2   given  by  (4  4)  and  S,  =  B/b,  ,   S2  =  B/b2  . 

Figure  7  shows  the  range  of  optimal  x2  as 
a  function  of  B  for  this  case.   It  is  immediately  evident 
that   x2  =  BQ/b,r,   is  optimal  regardless  of  the  value  of  B. 
(2)   B2  =  b2S  2  ,   B-j^  =  b1S  1   and   B1  <  B  <  B2  . 

By  similar  arguments  to  those  given  above  in 

*  *  it 

section  (1)  we  get   S,  =  S  ,   and   S2  =  B/b2  .   The  x2  range 
from  (33)  reduces  to 


43 


Xp4 


T  - 


Vbiri 

Bl/blrlt 


KdUlUulttllitttititleilw1 


\IUUIUU(ULU((IUU((UI 


litUUUftt 


FIGURE  7:     The  range  of  optimal   x?  as  a  function  of  B 
for  the  case  where     B     <  B~     and     B-.   <  Bp   . 


44 


S.  i  b,  S  , 

rl     Z  D2r2 


when   B,  <  B  <  B2  .   When   B  _<  B,  ,   S,  =  B/b,  , 

*  * 

S2  =  B/b2   and  x2  can  be  obtained  from  (44) .   Figure  7  also 

illustrates  the  range  of  optimal  x2  values  for  this  case. 

(3)   B,  =  b,S  -,  =  b0S  0  =  B0 
1     1  ul     2  u2     2 

In  this  case,  the  interval   B,  <  B  <  B2 

*  * 

does  not  exist  and  therefore  we  have   S,  =  S,,n  .   Sn  =  S  0 

1     ui  '    2     u2 

and  the  x2  interval  is  given  by  (37)  when   B,  =  B2  <  B  <  B^ 

■k  -k 

When   B  <_   B,   then   S,  =  B/b,  ,      S2  =  B/b2   and  the  x2 
interval  is  given  by  (44) .   Figure  8  shows  the  optimal  range 
of  x2  as  a  function  of  B. 

Having  now  examined  all  possible  cases  it  can  be 
seen  that  there  are  three  major  possibilities: 

First:   the  problem  is  unconstrained  if   B  j>  B-.   and 

partially  unconstrained  (the  range  of  x2  is  limited)' 

■k  * 

if  max  [B2,B  ]  <  B  <  B-  .   Solutions  for  S,  and  S2 

are  the  economic  order  levels  S  ,  and  S  0 . 

ul      u2 

Second:  If  the  storage  constraint  is  between  min[B,,BQ]  and 
max[B0,B  ]  then  the  solutions  depend  on  the  relation- 
ship  of  the  parameters  to  each  other  and  to  B.   They 
can  be  classified  according  to  two  ranges  for  B: 
namely,   BQ  <  B  <  B    and  B,  <  B  <  B2  . 
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FIGURE  8:     The  range  of  optimal   x2  as  a  function  of  B 
for  the  case  where     Bv  <  B9     and     B-,  ■  B«   . 
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•k  * 

Solutions  for  S,  and  S~  are: 


*      A         *      A 
1.    Sj^   =  S±      ;     S2   =   S2 


2.    Sx   =   B/b1   ;    S2   =   BQ/b2 


3.    s!   =   Bo/bl   ;    S2   =   B/b2 


4.    Sl   =   Sul   ;    S2   =   B/b2 


5.    S*   =   B/b1   ;    S*   =   Su2 


The  first  three  solutions  correspond  to  the   Bn  <  B  <  B 

* 
case  and  x2  can  only  have  one  value.   The  last  two  solutions 

correspond  to  the   B-,  <  B  <  B2   case  and  x2  can  take  on  any 

values  within  a  computed  interval. 

Third:   In  all  cases   where   B  £  min  [B,,B0 ]   the  space 

constraint  is  so  restrictive  that  only  one  item  can 

*  * 

be  stored  at  a  time.   Hence   S,  =  B/b,   and  S?  =  B/b2 , 

and  the  interval  for  x?  expands  as  B  decreases. 
The  details  of  the  optimal  solutions  of  the  two-item  case  are 
summarized  by  the  flowchart  presented  in  Figure  9. 

C.   THE  n-ITEM  CASE 

The  n-item  case  can  be  easily  formulated  using  the  same 
approach  as  for  the  two-item  case.   The  general  n  constraint 
equations  can  be  written  for  the  times  at  which  the  items 
arrive  using  the  max  functions. 


4  7 


Calculate 

s  1 

=  equation 

(28) 

Su2 

=  equation 

(29) 

B3 

=  equation 

(35) 

Bo  = 

Bx  = 
B,   = 


Calculate: 
equation  (36) 

equation   (34) 

equation   (30) 

equation  (45) 


FIGURE  9:     Flowchart  for  determining  S-j ,  S~  and  x~   for  the 
two-item  deterministic  order  level   model. 
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yes 


S1  =  B/b] 
S*  =  B/b2 
x2  in  range  (44) 


yes 


* 

s,  =  s  , 

ul 


S2  =  B/b2 

* 
x^  in  range  (47) 


no 


51  =  B/b1 
* 

52  =  Su2 
* 

x2  in  range  (46) 


FIGURE  9:  (continuation) 
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Calculate: 
S,   =  equation  (42) 
Sp  =  equation  (43) 


S,   =  B/b. 


1 

S*  = 
^2 

x*  = 


BQ/b2 
Vblrl 


yes 


yes 


S1  =  B0/b1 
S*  =  C/b2 
x*  =  BQ/b1r1 


* 
x2  =B0/blrl 


FIGURE  9:    (continuation) 
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The  principle  difficulty  involved  in  trying  to  solve  this 
problem  is  the  number  of  possible  combinations  which  must  be 
considered  because  of  the  various  max  functions.   Because  of 
its  complexities  the  n-item  case  will  not  be  addressed  in 
this  thesis. 

D.   EXAMPLES 

1.   Example  Number  1 

A  warehouse  stores  two  items  with  the  following 
parameters: 


Item  Number  1: 

r1   =  200  units  per  year; 

3 
b,   =    5  ft  /unit; 

c, ,  =  $  2  per  unit-year; 

c2,  =  $30  per  unit-year; 


Item  Number  2 : 

r2   =  250  units  per  year; 

b2   =    3  ft  /unit; 

c,2  =  $  1  per  unit-year; 

c22  =  $2  5  per  unit-year. 


Total  warehouse  space  available  is  600  ft  ;  the  scheduling 
period  T  is  one  year.   How  many  units  of  each  item  should  be 
ordered  and  when  should  the  orders  arrive? 

Following  the  steps  of  Figure  9,  we  first  calculate 

S  n  ,  S  0  and  B.,. 
ul'   u2       3 


S  ,  =  187.5 
ul 


S  0  =  240.4 
u2 


B3  =  1658  . 
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Since   B  =  600  <  B^  ,   we  must  calculate  BQ,  B  ,  B,  and  B2. 


Bn    =  429  ;   B   =  1229  ;   B,  =  721  ;   B~  =  937  . 
0  x  1  2 


We  see  that  B   >  B2   and  BQ  <  B   and  we  will  have  to 
calculate  S,  and  S2 . 


S1   =    106.4       ;       S2    =    165.6    . 


A  R  ~  B 

Now,    because      S,    <   r~~  =    120      and      S0    <   r—  =   200    ,      the 

1   bl  «  I 

optimal  values  for  S,  and  S2  are  S,  and  S2 .   The  corresponding 

x2  value  is  given  by  equation  (38) .   In  summary,  the  optimal 

solutions  are: 


*  * 

S,  =  106.4  units   ;    S2  =  165.6  units  ; 


x9  =  0.429  years 

Total  costs  TC,  using  equation  (25) ,  are  $  1125  per  year. 

2 .   Example  Number  2 

The  available  warehouse  space  is  the  same  as  in 
Example  number  1.   However,  the  parameters  of  the  items  are: 


52 


Item  Number  1 : 


Item  Number  2 : 


r,   =  200  units  per  year; 

3 

b,   =    5  ft  /unit; 

c,,  =  $  5  per  unit-year; 
c21  =  $100  per  unit-year; 


r~   =  250  units  per  year; 


12 


'22 


=    3  ft  /unit; 

=  $  1  per  unit-year; 

=  $  3  per  unit-year. 


Calculating  S  ,,  S  ~,  an^  B-, 


Sul  =  190  ;   Su2  =  187.5  ;   B3  =  1512  . 


Since  B  =  600  is  less  than  B~  calculations  of  B.,  B  ,  B, , 


and  B~  are  required. 


BQ  =  429  ;   Bx  =  1083  ;  B±   =    562.5  ;   B2  =  950  . 


Because   B   >  B2   and  BQ  <  B   we  have  to  calculate   S,   and 


S2  • 


S1  =  182.5   ;    S2  =  38.8  . 


B 


-  =  120  .   Sn  is  not  a  feasible  value  for  S, 
1    b,  1  1 


Because   S,  > 

The  optimal  feasible  solutions  are: 


S,  =  B/b,  =  120  units  ;   S2  =  BQ/b2  =  143  units  ; 


x2  =  B0/b,r,  =  0.429  years. 


Total  costs,  using  equation  (25),  are  $  1890  per  year. 
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E.   COMMENTS 

The  major  difference  between  the  model  examined  by 
Naddor  and  the  model  just  developed  is  that  the  different 
items  in  Naddor ' s  model  all  arrive  at  the  same  time,  while 
the  model  being  considered  in  this  paper  permits  time  phasing 
of  the  arrival  time  for  each  item.   This  time  phasing  allows 
a  higher  average  inventory  level  for  each  item  and,  conse- 
quently, a  reduced  average  shortage  level  and  a  lower  total 
average  annual  cost  for  the  system.   For  the  data  of  Example 
number  one  the  optimal  order  levels  using  Naddor ' s  model  are 
S,  =  51.3  units  per  year  and  S~  =  114.5  units  per  year, 
yielding  total  average  annual  costs  of  $2600  (using  Equation 
(25)).   Time  phasing  therefore  "saves"  $1475  per  year. 

In  addition  to  Naddor' s  model  being  more  expensive,  the 
solutions  are  not  always  feasible.   Naddor 's  analysis 
implies  that  a  negative  optimal  order  level  "creates" 
additional  available  space.   In  the  model  of  this  paper  this 
problem  has  been  carefully  considered  in  the  max  functions. 
Naddor  does  not  address  this  problem  at  alii   Using  Naddor' s 
results  to  solve  Example  number  2  we  get  S-,  =  177  units  and 
S2  =  -94  units.   Clearly,  b,S,  >  B  and  hence  S,  =  177  units 
is  infeasible. 
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V.   SUMMARY  AND  RECOMMENDATIONS 

A.   SUMMARY 

Chapter  III  provides  a  minimum  cost  (ordering  plus 
holding)  solution  to  a  multi-item  lot  size  inventory  problem 
subject  to  a  constraint  on  the  maximum  total  space  available 
for  storing  the  inventory.   The  variables  which  were  opti- 
mized were  the  scheduling  period  (time  between  orders  of  a 
single  item)  and  the  inter-order  times  between  the  different 
items.   The  assumptions  which  were  made  in  the  formulation 
of  this  model  were  a  deterministic  demand,  infinite  replen- 
ishment rate,  no  shortages  allowed  and  a  common  scheduling 
period  for  all  items. 

Chapter  IV  provides  a  minimum  cost  (holding  plus  back- 
order)  solution  for  a  two-item  system  under  periodic  review 
and  subject  to  the  same  space  constraint  as  in  Chapter  III. 
The  variables  which  were  optimized  were  the  order  levels 
for  the  two  items  and  the  inter-order  time  between  them. 
The  solutions  have  several  possible  forms  depending  upon 
the  interactions  of  the  constraint  and  the  parameters  of 
the  items.   The  assumptions  which  were  made  in  the  formula- 
tion of  the  problem  were  a  deterministic  demand,  infinite 
replenishment  rate  and  a  prescribed  interval  between  reviews 
which  is  identical  for  all  items. 
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B.   RECOMMENDATIONS 

Basic  to  the  models  presented  in  this  thesis  is  the 
assumption  of  deterministic  demand.   In  reality,  we  are 
faced  with  random  demand  patterns  and  therefore  it  is  impossi- 
ble to  determine  with  certainty  the  future  states  of  an  inven- 
tory system.   However,  it  has  been  shown  that  examination  of 
deterministic  inventory  models  can  provide  a  better  under- 
standing of  the  situation  being  modelled  and  may  suggest 
approximations  which  would  apply  in  the  case  of  random 
demand  [4],   This  thesis  provides  such  an  understanding  and 
one  of  the  next  steps  towards  improving  existing  severely 
constrained  multi-item  systems  is  to  obtain  demand  data  and 
classifications  of  items  as  either  needing  continuous 
reviewing  or  needing  only  a  periodic  review. 

The  analysis  of  the  lot  size  model  in  Chapter  III  was 
conducted  because  there  are  a  certain  number  of  items  which 
should  never  be  out  of  stock  in  a  BVP.   Further  research 
should  be  able  to  provide  a  listing  of  these  items  and  a 
probability  distribution  to  describe  their  underlying  demand. 
Because  these  items  should  be  continuously  reviewed  to  avoid 
stockouts  they  should  probably  be  stored  together  on  one  or 
several  trucks  to  facilitate  careful  monitoring. 

The  periodic  review  model  of  Chapter  IV  was  presented 
because  the  majority  of  the  items  which  have  to  be  stored  by 
a  BVP  are  not  of  such  vital  importance  that  they  cannot  be 
on  backorder  a  certain  amount  of  time.   This  model  emphasized 
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the  importance  of  knowing  the  values  of  time-dependent 
backorder  costs.   Such  costs  are  needed  and  military  essen- 
tiality coding  should  be  a  valuable  help  in  appraising  such 
costs  for  each  item,  thus  insuring  that  the  more  important 
items  have  higher  backorder  costs  and  will  therefore  be  less 
on  backorder  in  an  economic  inventory  system.   The  demand 
distribution  will  also  be  needed  before  the  model  results 
can  be  extended  to  a  random  demand  situation. 

Though  the  problems  of  solving  the  n-item  case  in 
Chapter  IV  seem  to  be  severe  for  large  n,  situations  involving 
large  n  values  probably  do  not  occur  in  a  BVP.   As  has  been 
already  mentioned,  a  BVP  consists  of  a  fleet  of  trucks. 
If  the  items  which  are  to  be  stored  on  one  truck  are  selected 
in  a  way  that  all  items  have  similar  parameters  —  low  cost 
and  high  demand  for  example  —  then  it  could  be  possible  to 
form  classes  of  items  and  treat  each  class  as  one  item.   Thus 
the  classes  per  truck  could  perhaps  be  reduced  to  two  or 
three.   Grouping  of  items  into  classes  of  similar  demand 
would  also  eliminate  the  need  to  determine  a  probability 
distribution  for  each  item;  the  demand  distribution  for  the 
class  would  suffice.   Of  course  the  notion  of  "similar  demand" 
is  rather  nebulous  and  would  require  some  further  defining 
before  grouping  could  be  made. 

While  demand  and  cost  data  are  being  gathered,  further 
analytical  work  should  proceed  using  the  ideas  behind  the 
approximate  random  demand  models  of  Chapter  4  and  5  of  Hadley 
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and  Whitin  [4].   In  particular  the  (Q,  r)  model  of  Chapter  4 
addresses  the  problem  of  optimal  ordering  when  backorder 
costs  are  high.   This  situation  corresponds  to  that  of  the 
BVP  for  items  which  should  "never  be  out  of  stock". 
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